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We employ the Unitary Correlation Operator Method (UCOM) to construct correlated, low-momentum matrix 
elements of realistic nucleon-nucleon interactions. The dominant short-range central and tensor correlations 
induced by the interaction are included explicitly by an unitary transformation. Using correlated momentum- 
space matrix elements of the Argonne V18 potential, we show that the unitary transformation eliminates the 
strong off-diagonal contributions caused by the short-range repulsion and the tensor interaction, and leaves 
a correlated interaction dominated by low-momentum contributions. We use correlated harmonic oscillator 
matrix elements as input for no-core shell model calculations for few-nucleon systems. Compared to the bare 
interaction, the convergence properties are dramatically improved. The bulk of the binding energy can already be 
obtained in very small model spaces or even with a single Slater determinant. Residual long-range correlations, 
not treated explicitly by the unitary transformation, can easily be described in model spaces of moderate size 
allowing for fast convergence. By varying the range of the tensor correlator we are able to map out the Tjon line 
and can in turn constrain the optimal correlator ranges. 
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I. INTRODUCTION 

One of the prime challenges in modern nuclear structure 
theory is the description of properties of nuclei across the 
whole nuclear chart based on realistic nucleon-nucleon in- 
teractions. Several modern nucleon-nucleon interactions that 
reproduce the experimental two-body data with high preci- 
sion are available, e.g., the Argonne V18 potential JlJ], the CD 
Bonn potential Q], or the Nijmegen potentials 01 . The use of 
these interactions for nuclear structure calculations in a strict 
ab initio fashion is restricted to light isotopes, where Green's 
Function Monte Carlo IH 01 or no-core shell model cal- 
culations 00,01 are computationally feasible. These virtu- 
ally exact solutions of the nuclear many-body problem show 
that realistic NN-potentials supplemented by a phenomeno- 
logical three-nucleon force are able to reproduce experimen- 
tal ground states and excitation spectra of light nuclei. Fur- 
thermore, recent developments in chiral perturbation theory 
provide schemes to construct two- and three-nucleon forces 
systematically JIoHTlll . 

A major obstacle for ab initio nuclear structure calcula- 
tions are the strong short-range correlations induced by re- 
alistic NN-interactions. There are two dominant components: 
(i) correlations induced by the short-range repulsive core in 
the central part of the potential and (it) correlations generated 
by the strong tensor interaction. It is well known that, in a 
shell-model language, the description of these correlations re- 
quires extremely large model-spaces — the repulsive core and 
the tensor interaction lead to sizable admixtures of high-lying 
shells. Simple many-body spaces, which remain tractable for 
large particle numbers, cannot describe these correlations. In 
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the extreme case, e.g. in a Hartree-Fock approach, the many- 
body state is restricted to a single Slater determinant which is 
not capable of representing these correlations by construction. 
Therefore, the use of a bare realistic NN-interaction in such a 
framework has to fail. 

There are several recent attempts to tackle this problem. 
One is the so-called Mowfc approach OB, which employs 
renormalization group techniques to reduce the bare realis- 
tic potential to a low-momentum interaction. Effectively, the 
high-momentum contributions, which are responsible for the 
admixture of high-lying states, are integrated out leaving an 
effective low-momentum interaction suitable for small model 
spaces. 

Another approach is the Unitary Correlation Operator 
Method (UCOM) Odd- 

Here the short-range central 
and tensor correlations are explicitly described by a state and 
basis-independent unitary transformation. Applying the uni- 
tary operator of the transformation to uncorrelated many -body 
states, e.g., the Slater determinant of a Hartree-Fock scheme, 
leads to a new correlated state which has the dominant short- 
range correlations built in. Alternatively, the correlation op- 
erator can be applied to the Hamiltonian, leading to a phase- 
shift equivalent correlated interaction Vucom which is well 
suited for small low-momentum model spaces. Hence it can 
be used as an universal input for a variety of many-body meth- 
ods. The operator form of this correlated interaction resulting 
for the Argonne V18 (AVI 8) potential has been used success- 
fully to perform nuclear structure calculations in the frame- 
work of Fermionic Molecular Dynamics fl^.fl7lfl^l . 

In this paper we are going to apply the Unitary Correla- 
tion Operator Method to derive correlated two-body matrix 
elements. They serve as convenient and universal input for 
a variety of many-body techniques, ranging from Hartree- 
Fock to shell-model. Following a summary of the formal- 
ism of the Unitary Correlation Operator Method in Sec. [H] 
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we derive explicit expressions for correlated matrix elements 
in Sec. [THI Optimal correlation functions for the AVI 8 po- 
tential are constructed in Sec. IIVI and the properties of the 
correlated momentum-space matrix elements are discussed in 
Sec. [V] Finally, in Sec. IVII we present results of no-core 
shell model calculations using correlated oscillator matrix ele- 
ments, which highlight the effect of the unitary transformation 
and the properties of the correlated interaction. 

II. THE UNITARY CORRELATION OPERATOR METHOD 

(UCOM) 

A. Unitary Correlation Operator 

The concept of the Unitary Correlation Operator Method 
can be summarized as follows: The dominant 
short-range central and tensor correlations are imprinted into a 
simple many -body state \^) through a state-independent uni- 
tary transformation 

|*)=C|tf). (1) 

The unitary correlation operator C describing this transfor- 
mation is given in an explicit operator form, independent of 
the particular representation or model space. The correlated 
many-body state explicitly contains the important short-range 
correlations generated by the interaction. Even if we start with 
a simple Slater determinant as uncorrelated state then the 
correlated state |W) cannot be represented by a single or a 
superposition of few Slater determinants anymore. 

When calculating expectation values or matrix elements of 
some operator A using correlated states 

(*| A |*') = (*| CUC |*') = (*| A |*') , (2) 

we can define a correlated operator through the similarity 
transformation 

A = C _1 AC = C^AC . (3) 

Due to the unitarity of C the notions of correlated states and 
correlated operators are equivalent and we may choose the 
form that is technically more advantageous. 

In the case of the nuclear many-body problem, the unitary 
correlation operator C has to account for short-range central 
and tensor correlations as outlined in Sec. U] It is convenient 
to disentangle these different types of correlations and define 
the correlation operator as a product of two unitary operators, 

C = CfiCy , (4) 

where Co describes short-range tensor correlations and C r 
central correlations. Each of these unitary operators is written 
as an exponential of a Hermitian two-body generator 



The construction of the generators g r and go, which encode 
the relevant physics of short-range interaction-induced corre- 
lations, is crucial. 

We start with the generator g r associated with the short- 
range central correlations induced by the repulsive core in the 
central part of the NN-interaction. At small relative distances, 
the two-body density is strongly suppressed as a result of the 
repulsive core. Pictorially, the core keeps the nucleons apart 
from each other so that they reside at larger distances outside 
the short-range repulsion fi^lUgl . These correlations can be 
imprinted into an uncorrelated many-body state by an unitary 
distance-dependent shift along the relative coordinate for each 
particle pair. Such radial shifts are generated by the projection 
of the relative momentum q = ^[Pi — p 2 ] onto the distance 
vector r = xi — X2 of two particles: 

qr = ^[f -q + q-f] • (6) 

The distance-dependence of the shift - large shifts at small 
distances within the core, small or no shifts outside the core - 
is described by a function sst(^) for each spin-isospin chan- 
nel. Their shape depends on the potential under consideration 
and contains all information on the short-range central corre- 
lations. The determination of the sgr(r) is discussed in detail 
in Sec. IIVI The full generator for the central correlations reads 

gr = ^2 o t SST ( r ) ^ + ^ s st (r)] n 5T , (7) 

S.T 

where Ust is the projection operator onto two-body spin S 
and isospin T. 

The correlations induced by the tensor part of the interac- 
tion are of a more complicated nature. They entangle the spins 
of the two nucleons with the direction of their relative distance 
vector r. Depending on the orientations of the spins, the nu- 
cleons are shifted perpendicular to the relative distance vector 
lfl~5l fl6ll - Such shifts are generated by the residue of the rel- 
ative momentum operator after subtracting the radial compo- 
nent 

r 1 

q n =q- -q r = ^(L x r-r x L) . (8) 

This "orbital momentum", embedded into a tensor operator 
which encodes the complicated entanglement between spatial 
and spin degrees of freedom, enters into the generator of the 
tensor correlations 

T 

using the general definition 

s 12 (a,b) =|[(<ri.a)(<r 2 -b) + (o-i-b)(<r 3 -a)] 
- i(cr 1 -cr2)(a • b + b • a) . 



Co = cxp [ - i gn >^ ' C r = exp |^ - i X g n ^ . (5) 

i<j i<j 



Note, the tensor operator s 12 (r, q n ) entering into the genera- 
tor go has the same structure as the standard tensor operator 
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S12 = Si2(f , f ) appearing in the bare potential except for the 
replacement of one of the relative coordinate vectors by the 
orbital momentum. Similar to the central correlators the func- 
tions t?r(r) describe the distance dependence of this angular 
shift for isospin T = and T = 1. Both, sst{t) and $t(j) 
have to be in accord with the potential under consideration. 

The crucial difference between the Unitary Correlation Op- 
erator Method and other schemes using similarity transforma- 
tions to construct an effective interaction, such as the Lee- 
Suzuki transformation lll9ll or the Unitary Model Operator 
Approach 12011 - is that our unitary correlation operator is given 
in an explicit operator form. This enables us to evaluate cor- 
related wave functions or correlated operators analytically as 
will be shown in the following. 



B. Correlated Wave Functions 

We consider the effect of the correlation operators on the 
component of a two-nucleon state that describes the relative 
motion. The center of mass part is not affected by the uni- 
tary correlators because they depend only on relative posi- 
tions and momenta. For the uncorrelated relative wave func- 
tion we assume LS'-coupled angular momentum eigenstates 
\4>{LS)J M TMt). For the sake of simplicity, the projection 
quantum numbers M and Mt are omitted in the following. 

The central correlator c r = exp(— ig r ) fHll affects only the 
radial part of the state and leaves the angular momentum and 
spin components unchanged. In coordinate representation it 
resembles a norm-conserving coordinate transformation []7" 



R-(r) 



(R-(r)\cf>) 
(r\4\<f>) = ?±^jR^f)(R + (r)\<t>) 



(11) 



where R+(r) and R-(r) are mutually inverse, R±[Rzp(r)] = 
r. These correlation functions are related to the function s(r) 
in the generator through the integral equation 



±1 



(12) 



To a certain approximation the following intuitive relation 
holds R±(r) « r ± s(r). For the sake of brevity we omit 
the spin and isospin indices of the correlation functions here 
and in the following. 

The action of the tensor correlator Co on LS'-coupled two- 
body states can be evaluated directly 1151 . The matrix ele- 
ments of the tensor operator Si2(r, q n ) for those states have 
only off-diagonal contributions 



{4>{J ± 1, l)JT\ sia(r, q ) \<i>{J T 1, l)JT) 
= ±3i y / J(J+l) . 



(13) 



Within a subspace of fixed J one can easily obtain the matrix 
exponential and thus the matrix elements of the full tensor 
correlator cq. On this basis we can construct explicit relations 
for the tensor correlated two-body states. 




FIG. 1 : (Color online) Application of the central and tensor correla- 
tors to a deuteron-like two-body wave function. Panels (a), (c), and 
(e) depict the uncorrelated, the central correlated, and the fully cor- 
related radial wave functions, resp. The panels (b) and (d) show the 
corresponding central and tensor correlation functions (see text). 



States with L = J are invariant under transformation with 
the tensor correlation operator 



co \4(JS) JT) = \<t>{JS)JT) 



(14) 



Only states with L = J ± 1 are susceptible to tensor correla- 
tions and transform like 



en \<f>(J ± 1, 1) JT) = cos 6j(t) \(f>{J ± 1, 1)JT) 
T sin^(r) |0(Jt1,1)JT) 

with r being the radial distance operator and 

6j{v)=3y/J{J + l) 0(r) . 



(15) 



(16) 



The tensor correlator admixes a state with AL = ±2 and 
changes the radial wave function of both components depend- 
ing on the tensor correlation function $(r). 

To illustrate the impact of the central and tensor correlation 
operators on a two-body state, we consider the example of 
the deuteron. Assume an uncorrelated state \<po(LS)JT) = 
|(/>o(01)10) which is a pure S*-wave state with the spin-isospin 
quantum numbers of the deuteron. The radial wave function 
(r\4>o) shall not contain short-range correlations induced by 
the repulsive core. Fig. [Qa) shows the uncorrelated L = 
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radial wave function. Applying the central correlator c r with 
the correlation function R + (r) depicted in panel (b) leads to a 
wave function which has the short-range central correlations, 
i.e. the hole at small interparticle distances, built in as de- 
picted in Fig. [flc). Subsequent application of the tensor cor- 
relator Co with the correlation function $(r) depicted in panel 
(d) generates the fully correlated wave function shown in Fig. 
[He). As a result of the tensor correlations, the wave function 
acquires a D-wave admixture whose radial structure depends 
crucially on the tensor correlation function. In order to repre- 
sent the long-range D-wave admixture, which is characteristic 
for the realistic deuteron wave function, a long-ranged ten- 
sor correlator is required as indicated by the dashed curve in 
panels (d) and (e). In the following sections we will argue that 
only the short-range and state-independent correlations should 
be described by the correlation operator. The solid curves in 
panels (d) and (e) correspond to an optimal short-range tensor 
correlator as will be constructed later-on (see Sec. IIVI optimal 
correlator for 1$ = 0.09 fm 3 ). 



C. Correlated Operators and Cluster Expansion 

The explicit formulation of correlated wave functions for 
the many-body problem becomes technically increasingly 
complicated and the equivalent notion of correlated operators 
proves more convenient. 

The similarity transformation Q of an arbitrary operator A 
leads to a correlated operator which contains irreducible con- 
tributions to all particle numbers. We can formulate a cluster 
expansion of the correlated operator, 



C t AC = AW+A[ 2 ]+A[ 3 l + 



(17) 



where A^ denotes the irreducible n-body part iM . When 
starting with a fc-body operator, all irreducible contributions 
Aj™] with n < k vanish. Hence, the unitary transformation 
of a two-body operator — the NN-interaction for example — 
yields a correlated operator containing a two-body contribu- 
tion, a three-body term, etc. 

The significance of the higher order terms depends on the 
range of the central and tensor correlations 11141 fl5l Hal . If 
the range of the correlation functions is small compared to the 
mean interparticle distance, then three-body and higher-order 
terms of the cluster expansion are negligible. Discarding these 
higher-order contributions leads to the two-body approxima- 
tion 



A C2 = A [i] + A [: 



(18) 



In principle, the higher-order contributions to the cluster ex- 
pansion can be evaluated systematically PUl . However, for 
many-body calculations the inclusion of those terms is an ex- 
treme challenge. 

Therefore, we restrict ourselves to the two-body approxi- 
mation and choose the correlation functions such that its ap- 
plicability is ensured. As discussed in detail in Sec. IVII we 
can use exact solutions of the many-body problem, e.g. in 
the no-core shell model framework, to estimate the size of the 
omitted higher-order contributions. 



D. Correlated Hamiltonian - Central Correlations 

In two-body approximation the unitary transformation of 
any relevant operator with the central correlation operator can 
be evaluated analytically. 

As a comprehensive example we consider a Hamiltonian 
consisting of kinetic energy and a realistic NN-interaction. 
For convenience we assume a generic operator form of the 
interaction 

V = E llM*)Op + °pV P (*)] (19) 



with 



O p = {1, (eri-er 2 ), q 2 , ql(cri -a 2 ), L 2 , L 2 (eri • er 2 ) 
(L-S), s 12 (f,f), s 12 (L,L)}®{1, (n-7-2)} 



(20) 



In order to accommodate momentum dependent terms, as they 
appear, e.g., in the Nijmegen or Bonn A/B potentials l22ll . 
we have chosen an explicitly symmetrized form. Notice, that 
any quadratic momentum dependence can be expressed by the 
q 2 and L 2 terms contained in d!9i . For simplicity, charge de- 
pendent terms are not explicitly discussed here. Nevertheless, 
they will be included in Sec. IVII 

For the formulation of the correlated Hamiltonian in two- 
body approximation, it is sufficient to consider the Hamilto- 
nian for a two-nucleon system given by 



h = T + v = t cm + t r + tj 2 + v , 



(21) 



where we have decomposed the kinetic energy operator T into 
a center of mass contribution t cm and a relative contribution 
which in turn is written as a sum of a radial and an angular 
part 



1 2 



1 L 

tfi — j 

2/1 x A 



(22) 



Applying the central correlator c r in two-body space leads to a 
correlated Hamiltonian consisting of the bare kinetic energy T 
and two-body contributions for the correlated radial and angu- 

— |"2l ~[21 

lar kinetic energy, t r and t n , resp., as well as the correlated 
two-body interaction vl 2 l 



zl h< 



T + t W + t 



M _u7P] , ~[2] 



(23) 



The explicit operator form of the correlated terms can be 
derived from a few basic identities. The similarity transforma- 
tion for the relative distance operator r results in the operator- 
valued function i?+(r) 



The unitarity c£ 
transforms as 



4 rc r = i?+(r) • (24) 
c" 1 implies that an arbitrary function of r 



4/(r)c r = /(ctrcv) = /GR + (r)) 



(25) 
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The interpretation of the unitary transformation in terms of 
a norm-conserving coordinate transformation r i— > R + (r) is 
evident. For the radial momentum operator q r one finds the 
following correlated form ltl4ll 

1 1 



C r C[ r c r 



(26) 



With this, we obtain the following expression for the square 
of the radial momentum, which enters into the radial part of 
the kinetic energy 



c£ q r c r 



1 r 

2 



1 



1 



i?V(r)2 



+ w(i) (27) 



with an additional local term depending only on the correla- 
tion function R + (r) 

, , 7R'i(r) 2 R'"(r) 



AR' + {r) 4 2R' + {rf ' 

All other basic operators, such as L 2 , (L-S), s 12 etc. commute 
with the correlation operator c r and are therefore invariant. 

Based on these elementary relations we can explicitly con- 
struct the two-body contributions to the correlated kinetic en- 
ergy. For the radial part we obtain using \21\ 



t 



1 



1 



1 



-2 Mr (r)^ ^2/i r (r) 
with a distance-dependent effective mass term 
1 If 1 



1 



(29) 



(30) 



2/i r (r) 2^\R' + (rf 

The two-body contribution to the correlated angular part of the 
kinetic energy involves only the basic relation (I26> and reads 

■ 2 



r[2] 



4tfic r 



tn 



1 



2Mn(r) r 2 

with a distance-dependent angular effective mass term 



1 



(31) 



(32) 



2/in(r) 2fi\R + (r)z 

The momentum dependent terms of the NN-interaction dl9l 
transform in a similar manner like the kinetic energy. Using 
(Eg) and J27j we obtain 

= l/^+(r)) 
2 V i?' 



q 2 

2 ^r 



n 2«(i2+W) 
4r 



i?' (rf 



(33) 



+ v(it f (r))u;(r)-« / (ii f (r)) 



R' + {rf 



For all other terms of the NN-interaction ( 1191 the commutator 
relations [q r , O p ) ~ [r, O p ] = are fulfilled and the similarity 
transformation with the central correlator reduces to 



4u(r) O p c r = v(R+(r))O p 



(34) 



Many of the other relevant operators, e.g. the quadratic radius 
or transition operators, can be transformed just as easily. 



E. Correlated Hamiltonian - Tensor Correlations 

The transformation of the Hamiltonian with the tensor cor- 
relation operator cn is more involved. In general, it can be 
evaluated via the Baker-Campbell-Hausdorff expansion 



c^Aco 



A + i[gn,A] + -[gn,[gn,A]] 



(35) 



Evaluation of the iterated commutators in some cases results 
in a termination of the series expansion. A trivial case is the 
distance operator r which commutes with the tensor generator 
go and is thus invariant under the transformation 



Cn r c n = r 



(36) 



For the radial momentum operator q r , the expansion J35I ter- 
minates after the first order commutators and we obtain the 
simple expression 



Cq qr co = q r - tf'(r) Si 2 (r, q n ) 



(37) 



Likewise, we find for the tensor correlated quadratic radial 
momentum operator 



Cq q 2 . Co =q, 2 - [tf'(r) q r + q r tf'(r)] Si 2 (r, q n ) 
+ [tf'(r)s 12 (r,q )] 2 , 



(38) 



where Si 2 (r,qo) 2 = 9[S 2 + 3(L • S) + (L • S) 2 ]. For all 
other operators of the interaction dl9> , that involve angular de- 
grees of freedom, the Baker-Campbell-Hausdorff series does 
not terminate. Through the commutators additional tensor op- 
erators are generated. For example, the relevant first order 
commutators are 

[gn,s 12 (f , f )] = itf(r)[-24ni - 18 (L • S) + 3s 12 (f, f )] 
[go,(L.S)]=h9(r)[-s 12 (qo,qo)] 
[go,L 2 ] =ii?(r)[2Bia(q n ,q n )] 
[go,s i2 (L,L)] = h9(r)[7si 2 (qo,qo)] , 

(39) 

where 

si2(q n , q n ) = 2r 2 Si 2 (qo, On) + s i2(L, L) - \ 812(7, f ) ■ 

(40) 

The next order generates higher powers of the orbital angular 
momentum operator, e.g. an L 2 (L • S) term, in addition. The 
resulting accumulation of new operators enforces a truncation 
of the Baker-Campbell-Hausdorff expansion at some finite or- 
der for the operator representation 1 16]. The basis representa- 
tion introduced in Sec. IHIIdoes not require this approximation. 



F. Correlated Interaction - Vucom 

Subtracting the uncorrected kinetic energy operator from 
the central and tensor correlated Hamiltonian in two-body 
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space defines the correlated interaction vucom in two-body 
approximation: 

vucom = c]x n hcnc r - T . (41) 

If we start from a realistic interaction which is given in an 
operator representation, e.g. the AVI 8 potential, then the cor- 
related interaction also has a closed operator representation 

vucom = r [v p (r)O p + O p v p (t)] , (42) 
p 

where 

P = {1, (oi-ofc), q 2 ., ql((T 1 -a 2 ), L 2 , L 2 (cti-ct 2 ), 
(L-S), s 12 (f,f), s 12 (L,L), 

si2(qn>qa)i qrSi 2 (r,q n ), L 2 (L • S), 
L 2 si 2 (q n ,q n ),...}® {1, (ri-r 2 )}. 

The dots indicate that higher-order contributions of the Baker- 
Campbell-Hausdorff expansion for the tensor transformation 
have been omitted. The terms shown above result from a 
truncation of the expansion J35i after the third order in go. 
For most applications the inclusion of these terms is sufficient 

m. 

The existence of an operator representation of vucom is 
essential for many-body models which are not based on a 
simple oscillator or plane-wave basis. One example is the 
Fermionic Molecular Dynamics model I23I l24ll which uses 
a non-orthogonal Gaussian basis and does not easily allow for 
a partial wave decomposition of the relative two-body states. 
Nevertheless, it is possible to evaluate the two-body matrix el- 
ements of vucom analytically, which facilitates efficient com- 
putations with this extremely versatile basis fl6l[T7lfT8Tl . 

As we have emphasized already, the operators of all ob- 
servables have to be transformed consistently. The uni- 
tary transformation of observables like quadratic radii, den- 
sities, momentum distributions, or transition matrix elements 
is straightforward given the toolbox acquired for the trans- 
formation of the Hamiltonian. The Unitary Correlation Op- 
erator Method owes this simplicity to the explicit state and 
representation-independent form of the correlation operators. 
In contrast, in many other approaches for the construction of 
an effective interaction, e.g. the Lee-Suzuki transformation 
0,0, Q]3] or the V[ ow fc renormalization group method fl2ll . the 
consistent derivation of effective quantities other than the in- 
teraction is a complicated and rarely addressed problem [25]. 

An important feature of vucom results from the finite range 
of the correlation functions sst(t) and i?t(?") entering into 
the generators. Since the correlation functions vanish at large 
distances — i.e., the correlation operator acts as a unit op- 
erator at large r — asymptotic properties of a two-body 
wave function are preserved. This implies that vucom is 
by construction phase-shift equivalent to the original NN- 
interaction. The unitary transformation can, therefore, be 
viewed as a way to construct an infinite manifold of realistic 
potentials, which all give identical phase-shifts. 



It is interesting to observe in which way the unitary trans- 
formation changes the operator of the interaction while pre- 
serving the phase-shifts. The central correlator reduces the 
short-range repulsion in the local part of the interaction and, 
at the same time, creates a non-local repulsion through the 
momentum-dependent terms. The tensor correlator removes 
some strength from the local tensor interaction and creates ad- 
ditional central contributions as well as new non-local tensor 
terms. Hence, the unitary transformation exploits the free- 
dom to redistribute strength between local and non-local parts 
of the potential without changing the phase-shifts. The non- 
local tensor terms establish an interesting connection to the 
CD Bonn potential, which among the realistic potentials is 
the only one including non-local tensor contributions l26ll . 



III. CORRELATED TWO-BODY MATRIX ELEMENTS 

Having introduced the basic formalism of the Unitary Cor- 
relation Operator Method, we can now derive two-body ma- 
trix elements of the correlated interaction vucom- We consider 
relative LS-coupled states of the form \n(LS)JM TM T ), 
with a generic radial quantum number n, relative orbital an- 
gular momentum L, spin S, total angular momentum J, and 
isospin T. The matrix elements of vucom thus read 

(n (LS) JMTMt I vucom \n'(L'S)JMTM T ) 

= (n(LS) JMTMt | c+c a hc n c r - T \n'(L'S)JMTM T ) • 

(44) 

The center of mass part of the two-body states is irrelevant for 
the unitary transformation, since the correlation operator only 
acts on the relative degrees of freedom of the two-body states. 
In the following derivations the projection quantum numbers 
M and Mt are omitted for simplicity. The formal framework 
discussed in the following is completely independent of the 
particular choice of basis, only the angular momentum struc- 
ture is relevant. 

The Unitary Correlation Operator Method offers different 
ways to compute these matrix elements. If we assume a NN- 
interaction of the form d!9> , then we can use the operator rep- 
resentation (142 \ of vucom and evaluate the matrix elements 
directly. This approach is computationally quite efficient. If 
one expands the radial dependencies of the individual opera- 
tor channels in a sum of Gauss functions, all radial integrals 
can be calculated analytically. The matrix elements of the 
additional tensor operators contained in vucom can be given 
in closed form as well. However, this direct approach relies 
on the truncation of the Baker-Campbell-Hausdorff expansion 
(1351 employed to evaluate the tensor correlation. 

In order to avoid this approximation for the tensor transfor- 
mation we apply the tensor correlators to the two-body states 
and make use of the exact expressions dl4> and (1151 . The 
central correlators will be applied to the operator as before, 
since a simple and exact expression for the central correlated 
Hamiltonian exists (cf. Sec. Ill Di . We formally interchange 
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the ordering of the correlations operators using the identity 

cJcJjhcnCr = (cjxjjc,-) 4hc r (4cnc r ) 



= c' Q cl hc r cn 



(45) 



with the "centrally correlated" tensor correlation operator 

csi = 4cnc r = exp[-ii?(i?+(r)) si 2 (r, q n )] . (46) 

The central correlator commutes with s 12 (r,q n ) and trans- 
forms therefore only i?(r), see Eq. i25\ . The tensor correlator 
cq acts on TS-coupled two-body states with L = J like the 
identity operator (cf. Sec. Ill Bi 



c n \n{JS)JT) = \n{JS)JT) . (47) 
For states with L = J ± 1 we have the simple relation 

c n |n(J T 1, 1) JT) = cos#j(r) |n(J =p 1, 1) JT) 



±sin6»j(r) |n(J± 1,1) JT) , 



where 



> J (r)=3V7(J+l)i?(il f (r)) 



(48) 



(49) 



Using these relations we can calculate the correlated two-body 
matrix elements exactly. 



A. Interactions in operator representation 

We first consider a bare potential given in the generic oper- 
ator representation ( 1191 and derive the correlated matrix el- 
ements for the local contributions of the form w(r)0 with 
[r, O] = [q r , O] = 0, which includes all operators of the set 
(I20i except for the q 2 terms. 

The matrix elements for L — L' — J are not affected by the 
tensor correlations, only the central correlators act according 
to \3A\ . In coordinate representation we obtain 

{n(JS)JT\ c];4«(r)Ocnc r \n'{JS)JT) = 

dru^j(r)u n ,j(r) v(r) ((JS)JT\ O \(JS)JT) , 

(50) 

where v(r) = v(R+(r)) is the transformed radial dependence 
of the potential. The w ri! L(r) = r0„ i(r) are the radial rel- 
ative wave functions of the oscillator basis or any other basis 
under consideration. For the diagonal matrix elements with 

L = L' = J =p 1 we get 

(n(J=Fl, ljJTlcJxJ^rpcncv \n'(J^fl, 1) JT) = 



Kj T i( r ) u n ',j T i(r) v(r) 

x [((J=f1,1)JT|0|(Jt1,1)JT) cos 2 6j{r) 

+ ((J±l, 1) JT| O |(J±1, 1) JT) sin 2 9j{r) 

± {(JTl, 1) JT| O (J±l, 1) JT) 2 cos 0j(r) sin 6j(r)] 

(51) 



with 9j(r) = 6j(R + (r)). Finally, the off-diagonal matrix 
elements for L = J =p 1 and L' = J ± 1 read 

(n(JTl,l)JT|4cJ 2 u(r)Oc n c r n'( J±l, 1) JT) = 
dr u n,j T i( r ) u n' ,.m(r) v(r) 
x [((Jt1,1)JT|0|(J±1,1)JT) cos 2 0/(r) 
- ((J±l, 1) JT| O |(J=Fl, 1) JT) sin 2 6j(r) 
T ((Jt1,1)JT|0|(Jt1,1)JT) cos0j(r)sin0/(r) 

± ((J±1,1)JT|0|(J±1,1)JT) sin 0/ (r) cos j (r)] . 

(52) 

Apart from the integration involving the radial wave functions, 
the matrix elements of the operators O in LS'-coupled angular 
momentum states are required. Only for the standard tensor 
operator O = Si 2 (-, -) the off-diagonal terms on the right 
hand side of Eqs. 15 li and j52t contribute. For all other oper- 
ators in J20i the off-diagonal matrix elements vanish, and the 
above equations simplify significantly. 

The effect of the tensor correlator is clearly visible in 
the structure of the correlated matrix elements i5l\ and 
J52l >. It admixes components with AT = ±2 to the states. 
Therefore, the correlated matrix element consists of a lin- 
ear combination of diagonal and off-diagonal matrix elements 
((TS')JT| O \ (L' S)JT). In this way even simple operators, 
like L or (L • S) acquire non-vanishing off-diagonal corre- 
lated matrix elements d52t . 

The momentum dependent terms of the potential \\9\ allow 
for an exact evaluation of the similarity transformation on the 
operator level. For the tensor correlated form of the operator 



(53) 



we obtain 



^v 9rCo = -^(r)q 2 + q 2 w(r)] + u(r)[i?'(r)si 2 (r, q )] 2 



[«(r)i?' (r)q r + q r tf' (r)v(r)] s X2 (r, q n ) 



(54) 



by using Eq. fl38t . Subsequent inclusion of the central correla- 
tions leads to the following expression for the diagonal matrix 
elements with L = L' = J in coordinate representation: 



n{JS)JT\ 44vg r cnc r \n'{JS)JT) 



v{t)w{t)-v'(t) ^ 2 

I [Kj( r ) u n',j( r ) + u 'n.A r ) u n'j(r)] ^ 



(55) 



where v'(r) = v'(R + (r)). As before, the tensor correlator 
does not affect these matrix elements and only the central cor- 
relations lead to a modification. For the diagonal matrix ele- 
ments with L = L' = J =p 1 the tensor terms contribute and 
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we obtain 

(n(J=Fl, l)JT\ 44v 9r c nCr \n'(JTh 1)JT) 



1 



v(r) w(r) + v(r) 9' j(r) z 



~, , R +(r) 



+ U 



n, J:fl 



(r) w n ^, /T i(r)] — 



w(r) 



(56) 



with 0' j(r) = 0'j(R + (r)). Likewise, we find 

(n(J=Fl, l)JT\ 4c4v gr c fi c r |n'(J±l, 1)JT) 
= ± /dr [<^, /Tl (r) <, iJ±1 (r) - <*, /Tl (r) M„',j±i(r)] 
^ v{r)9'j(r) 

(57) 
J =F 1 and 



for the off-diagonal matrix elements with L 
L' = J±1. 

The matrix elements for the correlated radial and angular 
kinetic energy can be constructed as special cases of the in- 
teraction matrix elements discussed above. By setting v(r) = 
l/(2/i r (r)) in Eqs. d55l to i5H we obtain the matrix elements 
for the effective mass part of the correlated radial kinetic en- 
ergy ( I29> . The matrix elements of the additional local poten- 
tial in d29l and the angular kinetic energy ( 13 U follow directly 
from Eqs. to <|52}. 



B. Interactions in partial-wave representation 

So far we have discussed interactions given in a closed op- 
erator representation of the form ( 1191 . However, many mod- 
ern interactions, e.g., the CD Bonn potential or recent chiral 
potentials, are defined using a non-local partial-wave repre- 
sentation. This makes it difficult to employ them within many- 
body models which do not allow for a partial-wave expansion 
of the two-body states JT6ll . Nevertheless, the calculation of 
central and tensor correlated matrix elements of the form (I44i 
is straightforward for those interactions. 

Consider a general non-local NN-potential in partial-wave 
representation. For simplicity we assume the potential given 
in a generic coordinate space representation 



dr r 2 / dr' r' 2 

^ L,L',S,J,T (58) 

r(LS)JT) v WS j T {r,r') (r'(L'S)JT\ , 



where M and Mj- have been omitted for brevity. Interactions 
given in momentum space can be easily transformed into this 
representation. 

For the construction of the correlated matrix elements we 
only need the expressions for correlated two-body states used 



in the previous section. For L = L' = J the tensor correla- 
tions are not active and we obtain 



n(JS)JT\ 44 v c Q c r \n'{JS)JT) 



drrU + {r) J dr' ' r"R+{r') u^j(r) u n ,,j{r') (59) 
x v JiJiS ,j,T(r,r') , 

where v LL >sjT(r,r') = v L l'sjt(R+(t), R+(r')). Due to 
the non-local character with respect to the relative coordinate, 
the metric factors TZ+(r) = \J R' + (r) R + (r) / r resulting from 
the transformation of the radial wave function remain. For the 
diagonal matrix elements L = V = J =p 1 of the non-local 
interaction (1581 we obtain 



(n(J=Fl, 1) JT\ 44 2 vcflc, K(J=Fl, l)JT) = 

drrU+{r) J dr' r TZ+(r') u^, /Tl (r) u n /,j^.(r') 
x [vj T i,j T isj,T(r, r') cos 9 j (r) cos 6j(r') 
+ ^J±a,J±i,i,J,T(r, r') sin 8j(r) sin#j(r') 
± vj^.,JH,i,J,T(r, r') cos 9,j{r) sin 9j{r') 
± wj±a,JTi.i,J.T(r, r') sin 9 'j(r) cos 9j(r')] 

(60) 

with 9j(r) = 9j(R+(r)). Finally, the off-diagonal matrix 
elements with L = J =p 1 and L' = J ± 1 read 

(n(J=Fl, 1) JT\ 4cl vc fi c r [n'(J±l, 1) JT) = 

drr72.+ (r) J dr' r' TZ+(r') u* njTl {r) u„',j±i(r') 

x [wjri,j±i4,./,T(r, r') cos 0j(r) cos 6j{r') 
- wj±a,JTi,i,J,T(r, r') sin 6j(r) sin0j(r') 
T «j T i,j T i4,XT(r, r') cos9j(r) sin9j(r') 
± «Jii,j±t,i,j,T(r,r') sin 6j{r) cos 0j(r')] . 

(61) 

For local interactions VLL'SJT{r 7 r') = v l ™ L , SJT (r)S(r — 
r') j (r'r) the metric factors TZ+(r) can be eliminated and the 
above equations simplify substantially. For technical reasons 
we use these expression also for the calculations with the 
AVI 8 potential discussed in the following. 



IV. OPTIMAL CORRELATION FUNCTIONS 

The Unitary Correlation Operator Method encapsulates the 
physics of short-range central and tensor correlations in the 
set of correlation functions s(r) and $(r). In this section, we 
discuss a scheme to determine these correlation functions for a 
given NN-potential. One important task is to isolate the short- 
range state-independent correlations from residual long-range 
correlations that should not be described by the unitary trans- 
formation but by the many-body state. 
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The most convenient procedure to fix the correlation func- 
tions is based on an energy minimization in the two-body sys- 
tem fl5ll . For each combination of spin S and isospin T we 
compute the correlated energy expectation value using a two- 
body trial state with the lowest possible orbital angular mo- 
mentum L. The uncorrected radial wave function should not 
contain any of the short-range correlations, i.e., it should re- 
semble the short-range behavior of a non-interacting system. 
In the following we will use a free zero-energy scattering so- 
lution (f>L (r) oc r L . One could just as well use harmonic oscil- 
lator wave functions, the difference in the resulting correlation 
functions is marginal. 

The correlation functions are represented by parametriza- 
tions with typically three variational parameters. The long- 
range part is generally well-described by a double-exponential 
decay with variable range. For the short-range behavior, sev- 
eral different parametrizations have been compared. For the 
AVI 8 potential, the following two parametrizations for the 
central correlation functions have proven most appropriate: 

R\(r) =r + a {r/pf exp[- exp(r//?)] , 

R + {r) = r + a [1 — exp(— r/7)] exp[— exp(r//3)] . 

Which of these parametrizations is best suited for a particular 
channel will be decided on the basis of the minimal energy 
alone. Note that rather than s(r), we directly parametrize the 
function R + (r), which enters into the expressions for corre- 
lated operators and matrix elements. For the tensor correlation 
functions the following parametrization is used 

$(r) = a [1 — exp(— r/7)] exp[— exp(r//3)] . (63) 

The 5 = channels are only affected by the central cor- 
relators. Their parameters are determined from the energy 
minimization within the lowest possible orbital angular mo- 
mentum state, i.e. L = 1 for T = and L = for T = 1, 
resp., 

-Boo = (<Ml0)10|4hc r |<Ml0)10) , 

i (64) 
E01 = (4>o (00)01 1 4 hc r |0 O (OO)O1) . 

For S = 0, T = 1 the minimization of E01 by variation of 
the parameters of the central correlation function is straight- 
forward. The resulting parameters are summarized in Table|l] 
For S = 0, T = the potential is purely repulsive and, there- 
fore, the energy minimization, for a negligible gain in energy, 
leads to central correlation functions of very long range. In 
order to avoid this pathology we employ a constraint on the 
strength of the correlation function defined through 

I R+ = J drr 2 (R + (r)-r) . (65) 

The value of this constraint on the central correlation function 
for the S = 0, T = channel is fixed to I R+ = O.lfm 4 in 
accord with the typical values in the other channels. 

For 5 = 1 the tensor correlations are active as well and 
we determine the parameters of the central and the tensor cor- 
relation functions simultaneously. For T = the energy is 
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TABLE I: Parameters of the central correlation functions R+(r) for 
the AVI 8 potential obtained from two-body energy minimization. 
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r [fm] 

FIG. 2: (Color online) Optimal central correlation functions R+ (r ) — 
r for the AVI 8 potential according to the parameters given in Tab. 
The curves correspond to the different spin-isospin channels: 
(S,T) = (0,1) ( ■), (1,0) ( ), (0,0) ( ), (1,1) 

defined by the matrix element with L = states 

Ew = (^(OlJlOlcJc^hcnCr |<M01)10) ■ (66) 

For T = 1 the lowest possible orbital angular momentum is 
L = 1. From angular momentum coupling we obtain 0,1, and 
2 as possible values for J. Therefore, we define the energy 
functional which is used in the minimization procedure as the 
sum over all three possibilities with relative weights given by 
2J+ 1 

En = i^i(ll)01|4hc r 1^(11)01) 

+ §<<Mll)ll|4hc r |<M11)11> (67) 
+ f(0 1 (ll)21|4cl J hc n c r |0i(ll)21) . 

As mentioned earlier, the long-range character of the ten- 
sor force leads to long-range tensor correlations. However, 
long-range tensor correlation functions are not desirable for 
several reasons: (;) The optimal long-range behavior would 
depend strongly on the nucleus under consideration. Hence, 
our goal of extracting the state-independent, universal corre- 
lations forbids long-range correlation functions, (ii) The two- 
body approximation would not be applicable for long-range 
correlators. (Hi) Effectively, higher order contributions of the 
cluster expansion lead to a screening of long-range tensor cor- 
relations between two nucleons throu gh th e presence of other 
nucleons within the correlation range 11611 . For these reasons, 
we constrain the range of the tensor correlation functions in 
our variational procedure. We use the following integral con- 
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TABLE II: Parameters of the tensor correlation functions $(r) for 
the AVI 8 potential with different values 1$ for the range constraint 
obtained from two-body energy minimization. 



straint on the "volume" of the tensor correlation functions 



h = drr z i?(r) . 



(68) 



The constrained energy minimization for the S = 1,T = 
and the S = 1,T = 1 channels with different values of the 
tensor correlation volume Ig leads to optimal parameters re- 
ported in Table im The optimal parameters for the central cor- 
relation functions change only marginally with the tensor con- 
straint. Therefore, we adopt a fixed set of parameters for the 
central correlators given in Table U 

The optimal central correlation functions for the AVI 8 po- 
tential are depicted in Fig. |2 In the even channels, the correla- 
tion functions decrease rapidly and vanish beyond r « 1.5 fm. 
The central correlators in the odd channels are weaker and of 
slightly longer range due to the influence of the centrifugal 
barrier. For the tensor correlation functions the constraints on 
the range are important. Fig. [5]shows the triplet-even (a) and 
triplet-odd (b) tensor correlation functions $(r) for different 
I#. Because the tensor interaction is significantly weaker for 
T = 1 than for T = 0, the tensor correlator for this chan- 
nel has a much smaller amplitude. The relevant values for the 
constraint 1$ are therefore smaller for the triple-odd channel. 

We stress that the range constraint for the tensor correla- 
tion functions has an important physical and conceptual back- 
ground. The Unitary Correlation Operator Method is used 
to describe state-independent short-range correlations only. 
Long-range correlations of any kind have to be described by 
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FIG. 3: (Color online) Optimal tensor correlation functions #(r) for 
different values of the range constraint 1$. (a) Correlation functions 
for T = with 1$ = 0.06, 0.09, and 0.12 fm 3 . (b) Correlation 
functions for T = 1 with I# = 0.01, 0.03, and 0.06 fm 3 . The 
arrows indicate the direction of increasing 



the model space employed in the solution of the many-body 
problem. By constraining the range of the tensor correlators 
we set the separation scale between short-range and long- 
range correlations. The optimal value for tensor constraints 
cannot be fixed in the two-body system alone, but requires 
input from few-nucleon systems. We will come back to this 
point in Sec. lVIl 



V. PROPERTIES OF CORRELATED MOMENTUM-SPACE 
MATRIX ELEMENTS 

A. Effect of the Correlators 

In order to illustrate the effect of the unitary transformation 
in more detail, we discuss relative momentum space matrix 
elements of the form (q(LS)JT\ vucom W(L'S)JT), where 
q is the relative two-body momentum. The calculation of cor- 
related momentum-space matrix elements is performed using 
the relations derived in Sec. [HI] with radial wave functions 
given by the spherical Bessel functions. 

First, we consider the full set of matrix elements in the 
(q, g')-plane for the lowest partial waves and compare the bare 
AVI 8 potential with the correlated interaction. The plots in 
Fig. ^depict the matrix elements for the 1 Sq and the 3 Si par- 
tial waves as well as for the mixed 3 Si — 3 D 1 channel (from 
top to bottom). The left-hand column corresponds to matrix 
elements of the bare AVI 8 potential, the center column to cor- 
related matrix elements using the central correlator only, and 
the right-hand column to the fully correlated matrix elements 
of vucom including central and tensor correlators. 
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bare central correlations central & tensor correlations 




FIG. 4: (Color online) Relative momentum-space matrix elements of the bare AVI 8 potential (left-hand column), the central correlated AVI 8 
potential (center column), and the fully correlated AVI 8 potential. The different rows correspond to different partial waves: So (top row), 
3 Si (middle row), and 3 Si — 3 Di (bottom row). The optimal tensor correlator for 1$ = 0.09 fm 3 is used. The red dots mark the plane of 
vanishing matrix elements. The momenta are given in units of [fm -1 ] and the matrix elements in [MeV]. 



The gross effect of the unitary transformation on the domi- 
nant S'-wave matrix elements depicted in the upper two rows 
of Fig. |4] is similar. In both cases the matrix elements of 
the bare interaction are predominantly repulsive except for a 
very small region at small momenta. The inclusion of the cen- 
tral correlator, accounting for correlations induced by the re- 
pulsive core of the interaction, causes a substantial change in 
the correlated matrix elements. In a region of low momenta 
q, q' < 2fm~ 1 the matrix elements become strongly attrac- 
tive. For larger momenta the magnitude of the matrix elements 
is reduced, outside a band along the diagonal the momentum 
space matrix elements practically vanish. Only within this 
band a moderate repulsion remains. The inclusion of the ten- 
sor correlator does not change the matrix elements in the spin- 
singlet channel. In the spin-triplet channel the addition of ten- 
sor correlations enhances the effect of the central correlations. 



The attractive matrix elements at low momenta are enhanced 
while the off-diagonal matrix elements are further suppressed. 

These matrix elements demonstrate the two major effects 
of the unitary correlators: (/) For the important L = partial 
waves, the low-momentum matrix elements become strongly 
attractive as a result of the proper treatment of the correlations 
induced by the repulsive core and the tensor part, (if) The 
off-diagonal matrix elements outside a band along the diago- 
nal are strongly suppressed. Hence the unitary transformation 
acts like a pre-diagonalization. 

The importance of tensor correlations is accentuated in the 
3 Si — 3 Di channel depicted in the bottom row of Fig. |4] 
For the bare interaction only the tensor part contributes to this 
channel and the matrix elements reveal strong off-diagonal 
contributions. In many-body calculations, e.g. in a shell- 
model framework, these off-diagonal matrix elements are re- 
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tion results from a transformation of longer-ranged compo- 
nents of the tensor interaction into operator channels which 
are accessible to uncorrelated S'-wave states (cf. Sec. Ill F> . 

The off-diagonal matrix elements in the 3 Si — 3 Pi chan- 
nel show a complementary behavior. The lower panel in Fig. 
[5] depicts the off-diagonal matrix elements as functions of q' 
for fixed q = 0. As mentioned earlier, the matrix elements 
far off the diagonal are strongly suppressed by the unitary 
transformation — they are associated with short-range tensor 
correlations. With increasing range of the tensor correlator, 
off-diagonal matrix elements at successively lower momenta 
are suppressed as well. Hence the band of non-vanishing ma- 
trix elements along the diagonal is narrowed with increasing 
correlator range. 



1 2 3 

q' [fin -1 ] 



FIG. 5: (Color online) Momentum-space matrix elements of the cor- 
related AVI 8 potential for different ranges of the tensor correlator 
(solid lines): 1$ = 0.06 fm 3 , 0.09 fm 3 , 0.12 fm 3 . The arrow indi- 
cates the direction of increasing correlator range. The dashed line 
represents the matrix elements of the bare potential. The upper panel 
shows diagonal matrix elements in the 3 Si channel as function of 
q = q . The lower panel depicts the off-diagonal 3 Si — 3 Z?i matrix 
elements for fixed q = as function of q' . 



sponsible for the admixture of high-lying basis states to the 
ground state contributing strongly to the binding energy l27ll . 
The effect of the central correlator on these matrix elements is 
marginal. The tensor correlator, however, causes a significant 
reduction of the off-diagonal matrix elements. Outside of a 
band along the diagonal the matrix elements vanish, i.e. the 
admixture of higher momenta or oscillator shells is suppressed 
significantly. 

The off-diagonal contributions from the tensor interaction 
are not fully suppressed by the tensor correlators — only 
the high-momentum components are eliminated. This cor- 
responds to the short-range part of the tensor correlations in 
coordinate space, which we constrained the tensor correlation 
operator to. 

The dependence of the momentum space matrix elements 
on the range of the tensor correlation functions is illustrated 
in Fig. |5] Note that only the tensor correlation functions are 
changed, the central correlators stay the same. Therefore, only 
the matrix elements in the spin-triplet channels change, and of 
those the S = 1, T = channels are affected most. The upper 
panel depicts the diagonal q = q' matrix elements for the 3 Si 
channel. With increasing correlator range as indicated by 
the arrow, the attraction at low momenta is enhanced. This 
can be easily understood in the picture of correlated states: 
Longer-ranged tensor correlators generate a longer-range D- 
wave admixture such that the tensor attraction of the bare po- 
tential can be exploited to a larger degree. In the picture of a 
correlated Hamiltonian, the increased low-momentum attrac- 



B. Comparison with V| 0W fc 

On the level of momentum-space matrix elements we 
can directly compare the correlated interaction vucom with 
the Vi ow k matrix elements resulting from a renormalization 
group decimation of the bare interaction mm. Both ap- 
proaches aim at the construction of a phase-shift equivalent 
low-momentum interaction, though their formal background 
is completely different. The V[owfc approach relies on a de- 
coupling of a low-momentum P-space, constrained by a mo- 
mentum cutoff A, and a complementary high-momentum Q- 
space via a similarity transformation. After a second transfor- 
mation in order to restore hermiticity, the momentum-space 
matrix elements within the P-space are obtained. Matrix el- 
ements between P and Q-space vanish by virtue of the de- 
coupling condition, Q-space matrix elements are discarded 
(hence violating unitarity) such that nonvanishing matrix el- 
ements exist only for momenta below the cutoff. In contrast 
to vucom the Mowfe approach is entirely formulated at the level 
of momentum-space matrix elements for the different partial 
waves. This entails that a general operator representation of 
the effective interaction is not directly accessible. 

Despite their formal differences the matrix elements of both 
methods show a remarkable agreement in the dominant par- 
tial waves. Fig. [6] compares the Mowfc matrix elements for a 
cutoff momentum A = 2.1 fm -1 with the momentum-space 
matrix elements of vucom obtained with the optimal correla- 
tors for I§ = 0.08 fm 3 . Up to momenta q « 1.5 fm -1 the 
matrix elements agree very well in most partial waves. For 
the dominant S'-wave channels the agreement extends right 
up to the cutoff momentum for Mowfc ■ Above the cutoff mo- 
mentum the V[owfc matrix elements are zero by construction. 
In contrast, the matrix elements of vucom continuously ex- 
tend to larger momenta. This reflects the different conceptual 
ideas: Whereas V\ ov ,k attempts a decimation of the interac- 
tion to low-momentum contributions below the cutoff scale, 
vucom uses a prediagonalization of the matrix-elements by a 
unitary transformation. 



13 




0.5 1 1.5 2 2.5 0.5 1 1.5 2 2.5 " 0.5 1 1.5 2 2.5 

q[fin [ ] q[{mr [ ] q [fiiT 1 ] 



FIG. 6: (Color online) Comparison of the diagonal momentum space matrix elements of the correlated interaction vucom (red lines) with the 
VJowfc matrix elements (blue dots) for the AVI 8 potential. The dashed line corresponds to the matrix elements of the bare potential. The optimal 
tensor correlators for 1$ = 0.08 fm 3 are used and the VJowfc momentum cutoff is A = 2.1 fm~ 1 . 



VI. FEW-BODY SHELL-MODEL CALCULATIONS 

The correlated interaction and the correlated matrix ele- 
ments can be used as input for all kinds of many-body calcu- 
lations. We have already discussed nuclear structure studies 
in the framework of Fermionic Molecular Dynamics 
which rely on the operator form of the correlated interaction. 
The correlated matrix elements serve as an input for shell- 
model or Hartree-Fock calculations. In a forthcoming pub- 
lication we will present nuclear structure calculations based 
on correlated realistic interactions in the framework of the 
Hartree-Fock and the Random Phase Approximation covering 
all mass regimes. 

In this section we discuss the application of the correlated 
matrix elements in no-core shell model calculations for 3 H 
and 4 He. These few-body systems provide important infor- 
mation on the correlated interaction beyond the two-body 
level. We use the no-core shell model code developed by Petr 
Navratil et al. ll28ll . It is formulated in a translationally invari- 
ant harmonic oscillator basis using Jacobi coordinates. Input 
for the shell-model diagonalization are the relative two-body 
matrix elements of the correlated AVI 8 potential, including 
charge dependent terms and Coulomb interaction. We stress 
that the Lee-Suzuki transformation usually employed in the 
no-core shell-model 0, l28ll is not used here. We only per- 
form a plain shell-model diagonalization. The task of trans- 
forming the bare interaction into an effective interaction suit- 
able for shell-model calculations in small model-spaces is per- 



formed by the unitary correlation operators. 

The effect of the unitary correlation operators in a no-core 
shell model calculation for the ground state of 4 He is illus- 
trated in Fig. Q For a given size of the model-space, charac- 
terized by the maximum relative oscillator quantum number 
A/"max = 2A^ max + L maj£ , the ground state energy is plotted as 
a function of the oscillator parameter hfl. The upper panel de- 
picts the shell-model result for the bare AVI 8 potential with- 
out any explicit treatment of correlations. All correlations in- 
duced by the interaction have to be described by the degrees of 
freedom of the model space alone. As expected, huge model 
spaces are required in order to adequately describe short-range 
correlations. Within the computational limits of 7V max < 16 
for 4 He, one is not able to achieve convergence with the bare 
interaction. The exact 4 He ground state energy for the AVI 8 
potential l29ll (marked by the horizontal line) is still somewhat 
lower than the result from the shell-model diagonalization for 

Mnax = 16. 

The convergence behavior changes dramatically once we 
use the correlated matrix elements instead of the bare ones. 
The lower panel in Fig. depicts the no-core shell model 
results for 4 He obtained with the correlated AVI 8 potential 
using the tensor correlator for I§ = 0.09 fm 3 in the domi- 
nant S = 1,T = channel. The S = 1,T = 1 channel 
is irrelevant for the nuclei considered in this section and the 
corresponding tensor correlation function is set to zero. The 
comparison with the calculation for the bare AVI 8 potential 
reveals three major effects of the unitary transformation: 
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FIG. 7: (Color online) Ground state energy of 4 He as function of 
the oscillator parameter hQ for different model model-space sizes 
A/"max = 0, 2, . . . , 16 as indicated by the labels. The upper panel 
shows results for the bare AVI 8 potential, the lower panel corre- 
sponds to the correlated potential vucom for h = 0.09 fm 3 . The 
horizontal lines re pres ent the exact binding energy for the bare po- 
tential taken from 1291. 



(i) The ground state energy for very small model spaces, 
e.g. jVmax = 0, for which the space consists of a single Slater 
determinant, is lowered dramatically. Evidently, the inclusion 
of the dominant short-range central and tensor correlations 
through the unitary transformation is sufficient to reproduce 
the bulk of the binding energy. 

(if) With increasing model-space size the energy is lowered 
by a moderate amount. The convergence is drastically im- 
proved. Fully converged results, featuring a flat energy curve 
over a significant range of oscillator parameters Ml can be ob- 
tained in spaces of moderate size. The energy gain compared 
to the results in small model-spaces can be attributed to resid- 
ual long-range correlations not described by the unitary corre- 
lator. In contrast to short-range correlations, these long-range 
correlations can be described quite easily in model-spaces of 
manageable size, hence the fast convergence. 

(Hi) The converged energy is generally below the exact 
ground state energy for the potential under consideration. This 
violation of the variational bound is solely due to the omission 
of the three- and four-body terms in the cluster expansion dl7> 
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FIG. 8: (Color online) Ground state energy of 3 H and 4 He as func- 
tion of the oscillator parameter hfl for the correlated AVI 8 potential 
(U = 0.09 fm 3 ) obtained in a no-core shell-model diagonalization. 
The different curves correspond to different model-space sizes A/'max 
as indicated by the labels. The horizontal lines represent the exact 
binding energies for the bare potential taken from 1 29]. 



of the correlated Hamiltonian. As a direct consequence of the 
unitarity of the transformation, the exact energy eigenvalues 
of the correlated Hamiltonian including all terms of the clus- 
ter expansion are identical to the exact energy eigenvalues of 
the bare Hamiltonian. Hence, the difference between the exact 
result using the bare interaction and the correlated interaction 
in two-body approximation equals the contribution of higher 
cluster orders and thus provides a quantitative measure for the 
quality of the two-body approximation. 

A more detailed view of the convergence behavior for the 
correlated AVI 8 potential is presented in Fig. [8] for 3 H and 
4 He. For both systems the three aforementioned effects can 
be observed. The convergence for 3 H is somewhat slower be- 
cause of the long-range structure of the wave function which 
cannot easily be described within the oscillator basis. 

The no-core shell model results directly reflect the ba- 
sic aims of the Unitary Correlation Operator Method. The 
short-range correlations are described explicitly by a state- 
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FIG. 9: (Color online) Binding energies of 4 He versus ! H. The blue 
discs show the results of exact Faddeev calculations obtained by A. 
Nogga et al. l29Tl using different modern NN-potentials. The green 
diamonds show results obtained by including simple three-body 
forces in addition to the realistic two-nucleon potentials [ 23] - The red 
triangles are the converged no-core shell model results for the corre- 
lated AV18 potential for different values I# = 0.03, . . . , 0.12 fm 3 of 
the range-constraint for the tensor correlator. 



independent unitary transformation such that the bulk of the 
binding energy can readily be obtained in very small model 
spaces. Residual system-dependent long-range correlations 
have to be described by the model space, which is easily pos- 
sible in the framework of the no-core shell model. 

Considering the omission of the higher-order terms in the 
cluster expansions, the no-core shell model results allow for 
a study of the range-dependence of these contributions. The 
difference between converged result and exact calculations 
with the bare potential, i.e., the size of the omitted three- and 
four-body terms, increases with increasing range of the ten- 
sor correlators. This dependence is summarized in Fig. |9j 
where we plot the ground state energies of 4 He and 3 H in the 
E( 3 H)-E( 4 He) plane. The results for different tensor correla- 
tors with range constraints 1$ =0.03, ...,0.12 fm 3 (cf. Tab. 
ITU are represented by the triangles. All those points fall onto a 
straight line. Moreover, this line coincides with the so-called 
Tj on-line 12^1 . which characterizes a correlation between the 
4 He and 3 H binding energies found for different realistic two- 
body potentials that reproduce the same phase shifts (blue 
discs in Fig. [9jl. It is not surprising that the correlated interac- 
tions follow the same trend, because the correlated interaction 
v ucom generates the same phase-shifts as the original poten- 
tial. For each set of correlators, the resulting correlated inter- 
action vijcom provides a new phase-shift equivalent realistic 
potential. The value of the range constraint 1$ can be used 
to map out the Tjon line. A similar behavior was observed 
for the Vi ow fc interaction as a function of the cutoff parameters 

m. 

In our calculation only the range of the tensor correlators 



in the spin-triplet channels is varied and of those only the 
S = 1, T = channel is relevant for 4 He and 3 H. Hence, the 
variation along the Tjon-line is related exclusively to the ten- 
sor correlations. As discussed in Sec. Ill El the unitary trans- 
formation of the Hamiltonian with the tensor correlators pro- 
duces additional momentum-dependent tensor terms in vjjcom 
as well as different central contributions. The non-local ten- 
sor contributions seem to play an important role regarding the 
Tjon-line: Increasing the strength of the non-local tensor con- 
tribution by increasing the strength of the tensor correlators 
shifts the binding energies towards the experimental values 
and away from the realistic potentials with purely local ten- 
sor contributions. This is in accord with the results for the 
CD Bonn potential, which is the only one among the high- 
precision NN-potentials including non-local tensor contribu- 
tions l2(Hl . 

By choosing an appropriate value for the range con- 
straint, one obtains a phase-shift equivalent two-body poten- 
tial, which produces binding energies for 4 He and 3 H close 
to the experimental point (cf. Fig. |9jl. This is the result of 
a subtle cancellation between different three-body contribu- 
tions. For exact calculations using the bare potentials one im- 
mediately finds that the two-body potential alone (blue discs 
in Fig. [5} does not generate enough binding. It has to be 
supplemented by a three-body force which produces a net at- 
traction. The green diamonds in Fig. |9]show results obtained 
by A. Nogga et al. |29] using various simple parametriza- 
tions of phenomenological three-nucleon forces supplement- 
ing the different realistic two-body potentials. For the calcu- 
lation with the correlated interactions we have not included 
any of the three-body contributions, i.e., neither the genuine 
three-body force nor the three-body contributions of the clus- 
ter expansion are taken into account. The proximity to the 
experimental point for, e.g. I# = 0.09 fm 3 , thus indicates, 
that the omitted three-body terms of the cluster expansion can 
be tuned such that they cancel the contributions of the genuine 
three-body force to a large extent. The influence of the gen- 
uine three-body force on the binding energies in these small 
systems can therefore be minimized by a proper choice of 
the correlators, i.e., by choosing the phase-shift equivalent 
two-body force which needs the weakest three-body force. 
This, however, does not mean that the three-body force can be 
avoided completely. One should keep in mind that the above 
observations refer to a single observable and to very small sys- 
tems only. 



VII. CONCLUSIONS 

The Unitary Correlation Operator Method provides a pow- 
erful and transparent tool to construct phase-shift equivalent 
low-momentum interactions by means of an explicit unitary 
transformation of a realistic NN-potential. The physics of 
short-range central and tensor correlations is encapsulated in 
the optimal correlators which are determined in the two-body 
system. For the long-ranged tensor component a separation 
of short-range state-independent correlations and long-range 
correlations is performed through an additional constraint on 
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the range of the correlators. Once the correlators are fixed, 
we can evaluate the unitary transformation of either states or 
operators directly. 

In the case of two-body matrix elements of the correlated 
Hamiltonian in an LS-coupled basis, it is convenient to map 
the unitary correlators onto the two-body angular momentum 
eigenstates. The resulting correlated matrix elements reveal 
some of the important features of the correlated interaction. 
The unitary transformation causes a pre-diagonalization of the 
Hamiltonian, i.e., large off-diagonal momentum-space matrix 
elements induced by the central core and the tensor interac- 
tion are eliminated and non-vanishing matrix elements remain 
solely in a band along the diagonal. 

Correlated matrix elements, e.g. with respect to a harmonic 
oscillator basis, serve as universal input for different many- 
body calculations l32ll . We have demonstrated the use of those 
matrix elements in the no-core shell model for A < 4. In 
comparison to a shell-model diagonalization with the bare in- 
teraction we observe a dramatic reduction of the ground state 
energy in very small model-spaces and a significant improve- 
ment of convergence with increasing size of the model space. 
Here the interplay between unitary correlator and model-space 
becomes evident: The unitary correlation operator describes 
the state-independent short-range correlations induced by the 
central and the tensor part of the interaction — this accounts 
for the bulk of the binding energy. State-dependent long-range 
correlations are described by the model-space — this leads to 
the moderate gain in binding energy with increasing model 



space size. The rapid convergence indicates that those long- 
range correlations, unlike the short-range correlations, can be 
quite easily treated in small and computationally accessible 
model-spaces. 

The no-core shell model calculations also provide a guide- 
line for the choice of the range-constraint for the tensor cor- 
relators. As function of the range of the tensor correlators 
a manifold of phase-shift equivalent potentials is generated 
which map out the Tjon line. The correlator range can be cho- 
sen such that the exact ground state energies for A < 4 are in 
good agreement with experiment. Thus the net impact of the 
residual three-nucleon force on the binding energies in these 
small systems can be minimized. 

The next step is to use the correlated matrix elements as in- 
put for nuclear structure calculations also for heavier isotopes. 
In a forthcoming publication we will discuss Hartree-Fock as 
well as RPA calculations across the whole nuclear chart using 
the same correlated realistic interactions. 



Acknowledgments 

We are very grateful to Petr Navratil for providing us with 
the no-core shell model code used in Sec. IVII We thank 
Achim Schwenk for the Vi ow fe matrix elements used in Fig. [6] 
This work is supported by the Deutsche Forschungsgemein- 
schaft through contract SFB 634. 



[1] R. B. Wiringa, V. G. J. Stoks, and R. Schiavilla, Phys. Rev. C 

51, 38 (1995). 
[2] R. Machleidt, Phys. Rev. C 63, 024001 (2001). 
[3] V. G. J. Stoks, R. A. M. Klomp, C. P. F. Terheggen, and J. J. 

de Swart, Phys. Rev. C 49, 2950 (1994). 
[4] S. C. Pieper, K. Varga, and R. B. Wiringa, Phys. Rev. C 66, 

044310(2002). 

[5] S. C. Pieper and R. B. Wiringa, Ann. Rev. Nucl. Part. Sci. 51, 
53 (2001). 

[6] S. C. Pieper, R. B. Wiringa, and J. Carlson, Phys. Rev. C 70, 
054325 (2004). 

[7] E. Caurier, P. Navratil, W. E. Ormand, and J. P. Vary, Phys. Rev. 

C66, 024314(2002). 
[8] P. Navratil and W. E. Ormand, Phys. Rev. Lett. 88, 152502 

(2002) . 

[9] P. Navratil, J. P. Vary, and B. R. Barrett, Phys. Rev. C 62, 
054311 (2000). 

[10] D. R. Entem and R. Machleidt, Phys. Rev C. 68, 041001(R) 

(2003) . 

[11] E. Epelbaum, A. Nogga, W. Glockle, H. Kamada, Ulf-G. 

MeiBner, and H. Witala, Phys. Rev. C 66, 064001 (2002). 
[12] S. K. Bogner, T. T. S. Kuo, and A. Schwenk, Phys. Rep. 386, 1 

(2003). 

[13] S. K. Bogner, T. T. S. Kuo, A. Schwenk, D. R. Entem, and 

R. Machleidt, Phys. Lett. B 576, 265 (2003). 
[14] H. Feldmeier, T. Neff, R. Roth, and J. Schnack, Nucl. Phys. 

A632, 61 (1998). 
[15] T. Neff and H. Feldmeier, Nucl. Phys. A713, 311 (2003). 
[16] R. Roth, T. Neff, H. Hergert, and H. Feldmeier, Nucl. Phys. A 



745, 3 (2004). 

[17] T. Neff and H. Feldmeier, Nucl. Phys. A 738, 357 (2004). 
[18] T. Neff, H. Feldmeier, and R. Roth, Nucl. Phys. A 752, 321 
(2005). 

[19] K. Suzuki and S. Y. Lee, Prog. Theo. Phys. 64, 2091 (1980). 
[20] S. Fujii, R. Okamoto, and K. Suzuki, Phys. Rev. C 69, 034328 
(2004). 

[21] R. Roth, Ph.D. thesis, Technische Universitat Darmstadt 

(2000) . 

[22] R. Machleidt, Adv. Nucl. Phys. 19, 189 (1989). 

[23] H. Feldmeier, Nucl. Phys. A515, 147 (1990). 

[24] H. Feldmeier and J. Schnack, Rev. Mod. Phys. 72, 655 (2000). 

[25] I. Stetcu, B. R. Barrett, P. Navratil, and J. P. Vary, Phys. Rev. C 

71, 044325 (2005). 
[26] R. Machleidt and I. Slaus, J. Phys. G: Nucl. Part. Phys. 27, R69 

(2001) . 

[27] J. P. Vary, P. U. Sauer, and C. W. Wong, Phys. Rev. C 7, 1776 
(1973). 

[28] P. Navratil, G. P. Kamuntavicius, and B. R. Barrett, Phys. Rev. 

C 61, 044001 (2000). 
[29] A. Nogga, H. Kamada, and W. Glockle, Phys. Rev. Lett. 85, 

944 (2000). 

[30] A. Nogga, S. K. Bogner, and A. Schwenk, Phys. Rev. C 70, 
061002(R) (2004). 

[31] Small (capital) letters are used for correlation operators in two- 
body (A-body) space. 

[32] An optimized code for the calculation of the correlated oscilla- 
tor matrix elements is available from the authors upon request. 



